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^ . Abstract 



The hadronic decays of quarkonium and the B meson inclusive decay into J/^ + X, 
if treated in the zero-binding approximation, suggest a value of the strong coupling 
constant much smaller than the value implied in the standard model by running from 
measurements at the Z. Thus, assuming the standard model is correct in the low 
energy region, there must be very substantial relativistic binding corrections to these 
processes. We discuss the wave function factor appearing in covariant treatments of 
quarkonium production and decay processes with special attention to the way in which 
the spin of the quarks carries over into the spin of the bound state. We find that 
Lorentz covariance requires that the bound state be a superposition of free quark 
and antiquark spinors notably different from the usual ones. We resolve a superficial 
apparent paradox suggesting that the relative momentum of the quarks should lie in 
a plane perpendicular to the spin quantization axis and we calculate the J/^ binding 
corrections to lepton pair decay and to the inclusive B decay. 
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1 Introduction 



Although perturbative QCD was originally seen as the source of the Zweig rule leading to 
the narrowness of heavy quarkonium , measurements of the strong coupling constant 
at the Z boson would imply, in the non-relativistic or zero-binding approximation (ZBA), 
hadronic decay rates significantly greater than those observed Thus, if the standard 
model is correct in the quarkonium region, corrections to these decay rates beyond the 
ZBA must greatly suppress these rates. A similar problem exists in the treatment of the 
inclusive B meson decay into J/\E' + X. After correcting an error in noted by 0, the 
color singlet model gives an as{mh) some 20% lower than suggested by extrapolating in the 
standard model from measurements at the Z boson (i.e., 0.154 ± 0.005 versus 0.20 ± 0.01). 
These low values of the strong coupling constant required in the ZBA together with low 
values obtained in other low energy measurements have stimulated interest in the light 
gluino hypothesis, but direct searches for such light supersymmetric particles have, up to 
now, all turned out negative. This increases the importance of a thorough investigation 
of quarkonium wave function effects beyond the ZBA. Quarkonium has a useful analog in 
positronium ||^ although the former case is complicated due to the non-Abelian nature of the 
binding force and due to non-perturbative confinement effects. There is continuing discussion 
of the much smaller binding effects in positronium to which we can compare our results. 
For definiteness we restrict our attention to the vector quarkonium states such as J/"^, 
T, and in particular to the J/\Ef. In a relativistic treatment the production amplitude of 
such a state with polarization e together with other particles from an arbitrary initial state 
takes the form 

A = e*^M,. (1.1) 

In the rest frame of the J/\&, the states of spin +1,-1, or zero along the z axis (quantization 
axis) are described by polarization vectors 

e'^(±)« = -^(0, ±1,^,0) 

e'^(O)^ = (0,0,0,1). (1.2) 

(Throughout this article we use the superscript R to indicate J/\I' rest frame values; for 
phase conventions see H], page 54.) 

The production of a J/\E' is, of course, equivalent to the production of a quark of mo- 
mentum p and spin A/2 together with an antiquark of momentum p and spin A/2 in a bound 
state wave function The J/"^ four-momentum, P, and the relative momentum, K, are 

P = p + p 

K = {p-p)/2 (1.3) 



with 



P^ = M' 

2 -2 2 



p = p = m (1.4) 



Thus 



= m"-— (1.5) 
4 

The momentum, K, is conjugate to the relative position of the quark and antiquark in the 
bound state wave function; we understand (|1.5|) , therefore, in the sense of an average over 
the wave function. The matrix element A in ( |1 . 1| ) is therefore related to a wave function 
factor, Fyj, which is linear in e* and is a matrix in the Dirac space of the heavy quarks: 

A = j ^^K^)^{K){F^UOp^. (1.6) 

Here a and /3 are spinor indices and F.^^ depends on the quark and antiquark momenta and 
spins while O carries the dependence on other variables in the problem. If the initial and 
final states contain heavy quarks only in the bound state, then the matrix element takes the 
form of a Dirac trace over F^ and other factors from O. For example, the decay amplitude 
for J/\E' into an electron positron pair is 

^ = / (0^^^^^^^^^ (Tr F^7m) u{p-)^^v{p^) ■ 2m, ■ Airaejq' (1.7) 
where is the virtual photon's squared four-momentum and 

Fyj = loFj^o- (1.8) 
We use dimensionless spinors normalized to 

uu = —vv = 1. (1.9) 

The J/"^ wave function is strongly affected by non-perturbative confinement effects and 
is therefore model dependent. In this article we ask what can be said assuming only that the 
wave function is spherically symmetric (an S wave), that the binding is small compared to 
the J/"^ mass, and that the bound state can be treated as a superposition of on-shell quark 
and antiquark. 

Since the production of a J/\E' is the production of a quark-antiquark in the appropriate 
wave function, we would expect that, if A = A, for suitable quark and antiquark spinors v 
and u, 

r2 2m — 

F^{X)aP = y -J^^'^iP^ \)aU{p, A)/3 (1.10) 

while in the case of A = —A we would have the spin projection zero quarkonium state 

FUOU = ^ • (1-11) 



It is understood that contains a unit matrix in the quark color space which, with the 
in eqs. |1.10| and |1 . 1 1| , puts the quark-antiquark pair into a color singlet state. The traces 
in eqs. |1.6| and include a color space trace. In the ZBA, the normalization factors are 
derived in and [1^. They carry over without change to the more general case. Here we 
have used the Fourier transform and normalization identities 



(2vr) 



X 



1. 



:i.i2i 



Assuming spherical symmetry and on-shell quarks, in addition to ( |1.5|) we have 

P P ' 



< K^K^ >= I g^^ 



M2 



;i.i3) 



Binding effects in any production or decay process can, therefore, be expressed as a power 
series in < >. 



2 The spin paradox 



We are free to choose any quantization axis along which to measure the spin of the J/\&. If 
we make the conventional choice of the z quantization axis, one might attempt to 



use the standard Bjorken-Drell []TT| spinors 

1 



u{p, A) 
v{p, A) 



2m{p^ + m) 
1 



with 



Xx 



2m{p^ + m 

(1 + AfTg) / 1 



(m + ^) 
(m — ^) 



Xx \ 

; 



icr2Xx 



In the case A = A we have 







iXX 0) = ip2(l + P3)(AfTi -i(T2)/4 



(2.1) 



(2.2) 



(2.3) 



2V2M 

where the R refers to the rest frame values of (|1.2| ) so that = M70. In ( p.3|) we have used 
the common decomposition of the four-by-four Dirac matrices into two two-by-two spaces 



7o 



P3, 



(2.4) 
(2.5) 



A 



Inserting ( |0| ) into ( |1.10| , |l.lID one would obtain 



{F^U = ^ iirn - ^)^^(A)(M + /'^)(m + ^)Jo)^^ (2.6) 



2M^3M(po + m)(Po + m) 

where /q is the unit matrix in the three- by-three quark color space. In general this naive use 
of ( |2.1| ) leads to a non-covariant Fy^, as can be seen by noting the explicit appearance of the 
J/\[' rest frame values of e and P. Problems with Lorentz covariance have been discussed 
by and attributed to non-perturbative effects. The non-covariance of (p.6|) is connected 
with a non-zero relative momentum K. The rest frame values of / are related to the values 
in a general frame where the J/\E' has four-momentum P'^ = (P°, P) by 



Only if i^' = does (|2.6| ) take the covariant form 

Although binding effects in QCD may be inherently non-perturbative, we see this problem 
not as one to be resolved by assuming non-perturbative effects, but as requiring consideration 
of what spinors can be covariant in a bound state problem. Going beyond the ZBA requires 
a more careful choice of the Dirac spinors u and v. Another symptom of the problem is 
associated with the way the spin of the quarks is carried over to the spin of the bound state, 
as we discuss shortly. 

The spin operator along a quantization axis n appropriate to act on a Dirac spinor of 
momentum p is 

Sn=^l,fi (2.9) 

provided p ■ fi = and = —1. Conventionally, one takes 

n^= (0,0,0,1), (2.10) 

i.e., pure spacelike in the z direction. However, this spin operator acting on the spinors of 
(|2.1|) gives the appropriate A/2 or A/2 only if 

P3=P3 = 0. (2.11) 

In the rest frame of the bound state this implies that 

7^3 = (2.12) 

but, if the J/\E' has unit spin along the z axis and there is no orbital angular momentum, then 
the constituent quark and antiquark must each have a unique spin 1/2 along this axis. This 
leads to the apparent paradox that the relative momentum, K, in the J/\l/ rest frame must 
lie in the plane perpendicular to the quantization axis, contrary to the requirement that an 
s-wave state must have a spherically symmetric distribution in the relative momentum. 



3 Recovering covariance; resolving the spin paradox 



A free fermion spinor can only depend on the particle momentum and the spin axis as in 
( |2.1| ). However, in order to recover covariance in the bound state and resolve the spin paradox 
one needs to employ not the usual Bjorken-Drell spinors ( |2.1| ) but instead 

u{p, A) 



^2m{M/2 + m) ^2M{Pq + M) 




viP,x)= , '"-^ , ''-^ V (3,i; 

^2m(M/2 + m) ^2M{Pq + M) \ ^^2X1 / 
With these spinors in place of (|2.1|) the wave function factor becomes 



) = m\ ——vip, Xjuip, X) = ; — . 3.2 

For the spin projection zero state, the central expression is the symmetrized combination as 
in ( |LTTD . 



Our result agrees with that found in a recent work ||12[ that has appeared since we began 
this research. While confirming their result, the current article emphasizes the covariance 
requirement of bound state fermion spinors depending not only on the fermion momentum 
but also on the bound state momentum. In addition we consider in detail how the spin 
structure of the constituents carries over to the spin of the bound state. A comparison with 
other earlier work is given below. 

In calculating a quarkonium production or decay process including first order corrections 
beyond the ZBA, the usual prescription is to expand the matrix element in the relative 
momentum and use the expectation values ( |1.5| , |1.13D as we will do in the next section for 
the quarkonium lepton pair decay rate and the B meson to J/\E' + X rate. The remainder 
of this section is devoted to justifying the above choice of spinors in the quarkonium bound 
state. 

One should note first that the quark and antiquark spin projections along any axis only 
need to sum to the quarkonium spin along that axis after performing the integration over 
relative momentum, K. For any fixed value of K there can be an orbital angular momen- 
tum that comes into the equation. Further, the usual spin states are defined by reference 
to axes in the particle rest frame, and axis directions are not Lorentz invariant. Bjorken 
and Drell |]ll| specified their spinors for nonzero fermion three-momentum in terms of the 



Lorentz transformation taking (m, 0,0,0) to {E,p), but they did not point out that such a 
transformation is not unique: any rotation in the rest frame can be applied first. Such a 
rotation is the difference between the spinors we propose in ( |3.1| ) and those of (^); such a 



rotation also is the difference between a single Lorentz boost and the same boost realized 
via two noncoUinear boosts. To construct the quark and antiquark spinors with the proper 
spin structure we first boost each from their differing rest frames into the rest frame of the 
J/\I', since the conventional use of eqs.( |1.2|) implies the choice of the z axis in the rest frame 



of the quarkonium for the defining quantization axis. 



Any spinor in the rest frame of a particle of mass m is boosted into the frame where this 
particle has four-momentum p'^ by the Lorentz transformation matrix 



T 

J- n 



m + 



2m(m + Pq) 



(3.3) 



Any covariant matrix, A, in Dirac space is boosted between these two frames by sandwiching 
the matrix in the particle rest frame between Tp on the left and 

m + 7o^ 



T =T 



^2m(m + Pq) 



(3.4) 



on the right. The spinors (|3.1| ) may be written by transforming the quark and antiquark 
spinors from their rest frame, where they have spin projections A/2 and A/2 respectively, first 
to the quarkonium rest frame where they have momenta p^ and p^, and then transforming 
to the general frame where the quarkonium has momentum P^. Thus 



u{p, A) 
v{p,\) - 



1 + A75^ 
-'p-'p^ 7, ""0 



1 + A75^ 
TpTpR vo. 



Here, 



Mo 



1 



V / 



(3.5) 



(3.6) 



Note that the spinors of ( PTT| ) are 



/ \ 



"^0 










1 



V -1 / 



A75^ 



2 

A75^ 



Mo 



-Vq. 



(3.7) 



(3.8) 
(3.9) 



The quark and antiquark momenta, p^ and p^, in the quarkonium rest frame are related 
to their values p and p in the general frame where the quarkonium has four-momentum 

^ (^pO^ p) by 

^ _ P.p_M 



p 



'R 



M 2 

P{po + M/2) 



p- 



Po + M 



(3.10) 



and equivalent equations for p^. 

The 7o's in the Tp and Tp of eqs. (|3^ ) commute with 75^ and give +1 and —1 respectively 
acting on uq and vq. We may then use 

Tp{m-^^)Tp = {m-^) (3.11) 



or 



Tp{m^i)^) = {m+f)Tp 

Tp{m-^^) = {m-^)Tp. (3.12) 



Again, the 70 in Tp acting on -Uq ot Vq gives ±1 respectively so the spinors of (^.5]) become 



identical to those of (|3.1j ). Combining v and u and using (|2.3| ) we have 



Fw = m\ ——vu = ——, -Tp == — Tpim+mln. 3.13 

V 3M (M/2 + 771)"^ Mv/T2M ^ > 

Then, using 

Tpf^{M + f ^)Tp = /^(M + f), (3.14) 



we have ( p.2| ). This result agrees with that of |T2|. One might ask, however, why the 



conventional spinors ( p.l|) lead to a non-covariant result while the spinors of ( p.l|) or 
lead to a satisfactory covariant result. The answer lies in the following observation. A 
quark spinor is given by a covariant factor times a rest frame spinor. The rest frame spinors 



combine to a covariant form in the meson rest frame (|2.3| ). From there one would have to 
boost to the general frame via Tp. But such a spinor would not satisfy the on-shell condition. 
The combination ( p. 3D sandwiched with T^r and T^r also provides a covariant form in the 
meson rest frame. This boosted by the Tp gives a covariant form in the general frame and 
also satisfies the Dirac on-shell conditions due to ( p.l2|) . 



To investigate how the quark and antiquark spins add to the quarkonium spin we consider 
the momentum dependent spin operators for quark and antiquark 

"^iiP) = -LpIpR 1 pRi P 

The Si{p) satisfy the SU{2) algebra for any p. In the quark rest frame we have 

^ = ^75^. = ^ (3.15) 
2 2m '^^ 2 ^ ' 

where Cj is a pure space-like unit vector along the i'th axis. Without forcing the relative 
momentum to lie in the x, y plane, Ss{p) and Ss{p) have the proper eigenvalues A/2 and A/2 
acting on the spinors of ( |3.1D . In addition we note that, in the quarkonium rest frame where 
Tp = l, 

Ssip) = Ssip^) = 1^751^^ (3.16) 



where 



m Pq + m 



= T,.^T,n = ^+ "1^ H (3-17) 



Pq, of course, is merely M/2. After averaging over K directions fl^ is in the ^ direction. For 
on-shell quarks fl^ is, after averaging, precisely ^. Thus this construction insures that in the 
quarkonium rest frame the quark and antiquark spins are taken with respect to the same n 
direction after averaging over relative momenta. It is this property that ensures covariance 
of the vu bound state and ensures that the spins of the quark and antiquark carry over 
properly to the spin of the quarkonium. 

In the general frame, where the quarkonium has four-momentum (P°, P), 



where 



Ssip) = 7^l5fi (3.18) 
2m 



= n^p - l^Q^^ (3.19) 
m 



and is the dimensionless four-vector 



mP'^ + Mp'' 

= Mw^y 

Tip is the four- vector that is n in the quarkonium rest frame. 

np = A(P)n = ( fi + — -i-^ I (3.21) 
where A(P) is the Lorentz transformation matrix that relates the quarkonium rest frame to 



the general frame. The four-dimensional spin axis n for the antiquark is given by eqs. |3.19 
and |3.20| replacing p by p. Only in the K = limit or after averaging over K are these axes 
equal. Note that p ■ n = p ■ n = P ■ np = as required by Lorentz invariance. Using these 
orthogonalities and the fact that 



it is easy to show that 



2 

where the u and v spinors here are those of (|3.1|) . Therefore 



(m + ^)75(^(M+ f) = (3.22) 



S3{p)u{p,\) = ^m(p,A) 

SsiPHpJ) = '^v{pj) (3.23) 



S3{p)F^ + F^Ssip) = ^^F^- (3.24) 



This makes manifest the addition of the quark and antiquark spin projections to compose 
that of the meson for arbitrary values of K. 

We turn now to a comparison of our result with that of other authors. We convert all 
results to those appropriate to production of a quarkonium state as opposed to a decay 
process and we ignore normalization factors. 

In one of the seminal early papers on quarkonium processes 0, there appears the ex- 
pression 

F^~(m-^)^(7° + l)(m + ^). (3.25) 

These authors did not attempt to treat the system beyond the ZBA and noted that the above 
result held in the bound state rest frame, but did not deal with the apparent non-covariance 
of the result. 



In [|T3|, the wave function factor was given as 

^ {m-f)f{m+ (3.26) 

This paper omitted the factor, (M + /'), which is important in going beyond the ZBA. The 
same form appears in [|I^ . 

In [^, the form 

F^^f{m+r) (3.27) 

is given without attempting an analysis beyond the ZBA. 

In 0, in addition to the form ( p. 26 ), there appears the non-covariant form 

^(7° + l). (3.28) 

Again, these authors did not deal with the apparent non-covariance, preferring to restrict 
their considerations to the meson rest frame. 

Repeating our analysis above, it is straightforward to write the wave function factor for 
a bound state of different quark species. For a quark of mass m and antiquark of mass m, 
the appropriate generalization is 

M^J12M{pq^ + m) (po^ + m) 
Here the quark and antiquark energies in the rest frame of the bound state are 

R M — rrf 



2 2M 

Po^ = IT - -^TV^- (3.30) 



2 2M 

In the case of non-equal masses, the 6 function in ( [L.6|) is modified to 



■ K — rn? ^ , ^ 

'^—1 A^' 



-I r\ 



4 Decay rates beyond the zero-binding approximation 



We would now like to use the bound state formalism to calculate the decay rates J/"^ ^ 
and B J/\l/ + X including the first corrections beyond the ZBA. In the decay of a J/^ as 
opposed to the production one uses the conjugate wave function factor ( |1.8|) . Using ( |1.10|) 



in the trace of ( |1.7| ) and averaging over the relative momenta using ( |1.5|jl.l3| ) we have 



Tr(F..,) = 1 - 3p„^,,). U yi2. J 1 + 3^ . (4.1) 



Thus, to lowest order in K"^, 



Using the Particle Data Group ||r3| charm quark mass m = (1.25 ± 0.1) GeV and J/^E' mass 
M = 3.097 GeV we see that the first order correction beyond the ZBA decreases the leptonic 
pair decay rate of the J/\l/ by a factor of 0.88 ± .04. This correction has been also found in 
[0. The effect goes in the wrong direction to explain the suppression of the hadronic J/\& 
decays relative to the leptonic pair decay. However, important contributions come from the 
binding corrections to the three gluon decay. Although the first order correction does reduce 
the three gluon decay rate, the authors of |jl2| note an apparent slow convergence of the 



expansion indicating that further study is needed. 

Turning to the decay B — » J/\l/ + X, the Dirac structure of the invariant amplitude is 



A ~ Mp.)YWPMp.) (-9.. + ^) (4.3) 



where Pl is the left handed chiral projection and 

k^=Pb,,-p,,=Pst^+Pf,- (4.4) 
It is easy to see that, to lowest order in and to first order in the decay rate is 

modified by the identical factor (1 + |^) that occurs in lepton pair decay. Thus, measured 
relative to lepton pair decay or, equivalently, using the leptonic rate as a measure of |\l/(0)p, 
the relativistic corrections do not change the inclusive J/\l/ decay rate of the B meson at 
zeroth order in the strong coupling constant. 

In summary, we have shown how the recovery of Lorentz covariance requires the use of 
quark spinors in the bound state that depend not only on the quark momentum but also 
on the bound state momentum. We have also noted how the relevant spin operators for 
quark and antiquark refer to momentum dependent axes that are only equal in an average 
sense. We have confirmed earlier recent results for the wave function factor and lepton-pair 
decay mode of quarkonium, and have extended the beyond-ZBA analyses to the inclusive 
J/\l/ decay of the B meson. 
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